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A DIGITAL HIGHER ORDER INTERPOLATION 
PATH CONTROLLER* 

By William Granville Batte 
Langley Research Center 

SUMMARY 

A two-axis digital higher order interpolation path controller for generating a smooth 
incremental path through discrete Cartesian input data is designed, constructed, and 
evaluated. 

The novel interpolation technique exhibits the following features: 

(1) Its digital implementation is simple relative to classical higher order interpo- 
lation schemes. 

(2) It is readily adaptable to incremental techniques and to the calculus of finite 
differences. 

(3) It accommodates unequal argument spacing. 

(4) It processes multivalued and closed contour functions. 

(5) It can accommodate raw or nonpreprocessed data. 

(6) It is readily expandable to multidimensions. 

Although the controller may be applied in many areas including control of machine 
tools, navigation, simulation, function generation, remote control, and automatic plotting, 
the experimental model is evaluated with a standard incremental x-y plotter as the out- 
put device and a punched paper tape reader as the input device. 

Laboratory tests were made on the implemented system and actual copies of its 
output are included. These results show that the approach is feasible and that the 
research objectives are met. 


INTRODUCTION 

The subject controller is a two-axis digital higher order interpolation device which 
generates a smooth incremental path through discrete, plane Cartesian input data. These 
data, which may fall into any one or more of the four quadrants, are sequentially ordered 

The information presented herein was included in a thesis submitted in partial 
fulfillment of the requirements for the degree of Doctor of Philosophy, Case Institute of 
Technology, Cleveland, Ohio, 1965. 



as received so that the generated path follows this ordering. Since no attempt is made to 
fit the data to some analytic function, the controller may be viewed as a device for auto- 
matically "applying the draftman's french curve." 

Although the controller is developed for use in automatic plotting, it may be applied 
in many other areas including control of machine tools, navigation, simulation, function 
generation, and remote control. Also, the mathematical scheme may be used as an inter- 
polation technique for certain types of general purpose computing. 

The controller is the result of research directed toward satisfying a need in auto- 
matic digital plotting. Automatic x-y plotters have replaced manual methods in most 
facilities where large quantities of data are displayed graphically. Most of these plotters 
produce point plots in which a pen or symbol head marks the locations corresponding to 
the various discrete x-y input values. In addition, some plotters (e.g., the Beckman 210 
tape -to -plotter system) can also produce continuous curves in which the discrete input 
data are connected by straight-line segments. Not one plotter, however, is effective in 
"placing the french curve through the points" as was commonly done in the manual 
methods. Consequently the primary objective of this research is to evolve a simple 
controller capable of supplying the "french curve" even under conditions of (a) data with 
unequal argument spacing, (b) multivalued or closed contour functions, and (c) raw or 
nonpreprocessed data. Secondary objectives require the hardware to be relatively simple 
and fast enough to keep the output device operating at its inherent maximum speed. 

Numerous digital controllers are in existence but each fails in one or more ways 
to satisfy these objectives. For example, Mergler’s machine tool controller (ref. 1) 
requires rate information to be supplied as a part of the input data and, hence, does not 
satisfy condition (c). His proposed quadratic (second order) technique recognizes its 
inadequacies for certain "geometry" and, in addition, requires "equal intervals of the 
argument," each of which renders the scheme inadequate for the present application. 

Several second-order systems are in existence but they too do not satisfy condi- 
tion (c). For example, the Bendix Dynapath system (ref. 2) requires, for circular inter- 
polation, that the initial and final positions of the radius vectors be supplied as a portion 
of the input data. Similarly, the Fuji system (ref. 3) requires that the input data be 
expressed in a coordinate system whose origin is at the center of the arc. Even more 
serious, however, is that the present application requires at least a third-order interpo- 
lation as demonstrated in reference 4. 

Ninke (ref. 5) recently adapted the Newton-Gauss interpolation technique to a digital 
third-order controller. Using digital differential analyzer (DDA) methods, the system 
consists of three digital integrators cascaded so that the third derivative introduced at 
the first integrator is successively integrated three times, this integration yielding the 
path function at the output of the last integrator. The initial conditions of the integrators 
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and corrections thereto must be supplied as a portion of the input data; hence, condi- 
tion (c) is violated. Also the simplicity with which these quantities are computed is a 
result of the assumption by Ninke that the input data are spaced at equal increments of 
the arguments; hence, condition (a) is also violated. 

A similar scheme is outlined in reference 6 where again the approach is based on 
equal arguments. 

As with most digital systems the controller has its analog counterpart. In many 
situations, however, the accuracy requirement renders the analog approach ineffective. 
Consider, for example, even simple linear interpolation where it is required to connect 
two points with a straight line. This problem is approached with the conventional analog 
closed-loop plotter by integrating the inputs to the two axes with the same effective time 
constants. The solution sounds simple, but the word "effective" requires a considera- 
tion of the unbalanced inertias and frictions on the two axes, the unbalanced breakaway 
torques, and, of course, settings of gain, damping, and so forth. Furthermore, higher- 
order interpolation is even more complex! 

The study of existing controllers and classical interpolation techniques as applied 
to the performance requirements of the subject controller leads to the following broad 
observations which are presented in the nature of hypotheses upon which the controllers 
design is based: 

(1) It is not important to fit the input data to some analytic function (e.g., a third- 
degree polynomial or perhaps some french curve logarithmic function); instead, it is 
more important to utilize some hardware-oriented (though perhaps less mathematical) 
scheme such as a goal seeking system. 

(2) Since the entire path must be generated, as opposed to merely evaluating a few 
isolated function values, a finite difference technique should be utilized such that the 
function change is computed and used to update the function along the generated path. 

(3) Since the finite difference technique will no doubt result in many isolated arith- 
metic units, the internal coding should be natural binary with negative numbers repre- 
sented in two’s complement notations. This coding will simplify the designs of the units. 

(4) Since any closed contour path contains at least two infinite slopes, a technique 
which obviates the large slope problem (refs. 5 and 6) must be used. 

(5) For simplicity of the experimental hardware, an incremental, open-loop control 
system should be used. (Where applications require absolute or incremental closed-loop 
systems, the experimental design is easily adapted to them.) 

(6) Since the path must be generated from any data point through the next point, the 
direction of the latter point with respect to the present position of the generated path is an 
important variable and should probably appear explicitly in the basic interpolation scheme. 
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Other philosophies peculiar to various parts of the controller are found in the 
appropriate design sections. 


SYMBOLS AND NOTATIONS 


Arithmetic: 


Literals 


a,b,c,d,e 

AJ 

(AJ) 

(AJx),(AJy) 

CJ 

(CJ) 

d^,dx,. . . 

d£ 

ds 

F 

K 

k,k-l,k+l 

L 

P 

Q 

R 

RJ 

(RJx),(RJy) 


set of input data points 

accumulator J where J is any integer 

quantity in AJ 

quantity in x or y portion of AJ 
counter J where J is any integer 
quantity in CJ 
infinitesimal differentials 

instantaneous curvature of a trajectory 

multiplier associated with binary rate multiplier 

constant which relates two equivalent angles and i/'g 

present, previous, and next values of a discrete quantity 

length of (i.e., number of bit positions or stages in) accumulator, counter, 
or register 

point on generated path 

proportionality constant which determines trajectory 

carriage return character on tape 

register J where J is any alphabetical character 

quantity in x or y portion of RJ 

trajectory path 
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S[] sign of quantity within brackets 

SWJ switch J where J is any integer or alphabetical character 

x,y,z Cartesian coordinates 

a tangent angle of generated path at last data point 

j8 tangent angle of generated path at next data point 

y corrected error distance between grid corner and s-curve — grid- line 

intersection 

6^C4 c ) change in coarse portion of quantity in counter 4 

6i//,6x,. . . finite difference or change in i|/, x, ... 

A "space" character on tape separating x- and y-data 

Ai^,A0,. . . total change in \f/, (p, ... 

e error between (p and 9 C 

0 motion angle required for coincidence with s-curve at next grid crossing; 

subscript c or f denotes coarse or fine portion of 6, respectively 

9 C( p merging function which starts from 0 C and ends at <p 

p predicted error distance between grid corner and s-curve — grid- line 

intersection 

ip steering angle or instantaneous angle of tangent to theoretical trajectory 

xjy angle of present position with respect to next data point 

| CJ | absolute value of (C J) 

Prime with symbol indicates quantity is expressed in transformed coordinate 
system. 
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Logic: 


A,B,C,. . . Boolean variables 

Ak’Ak + i augend (minuend) bit in accumulator before and after arithmetic operation, 
respectively 

Bjj addend (subtrahend) bit input to accumulator 

A[AJ].,B[AJ] A and B inputs of AJ 

N accumulator (counter) is subtracting (counting backwards) 

N[AJ] operation of AJ is negative 

Ir'k’^k+l carry (borrow) bit in accumulator before and after addition (subtraction), 
respectively 

P[AJj propagation signal P of AJ 

SPGJ shift pulse generator J where J is any integer 

Sj,Rj,Tj set, reset, and trigger (complementary) inputs to flip-flop J where J 
is any alphabetical character 

S[] sign of quantity within brackets is negative 

X "don't care” condition 

Qfj flip-flop J is setting (i.e., is undergoing ”0" to "1" transition) or logical 

signal J is turning on (J is any alphabetical character) 

Q!j T defined as aj except signal is delayed by nr units and J is in terms 

of positive logic 

0 angle 0 is being computed; 0 implies that angle j// is being computed 

0j flip-flop J is resetting (i.e., is undergoing "1" to "0" transition) or logical 

signal J is turning off (J is any alphabetical character) 


6 



I 


1 6x | , J 6y | 1 increment of motion is occurring in x- or y-direction (without regard 

to its sense) 

r unit of delay 

A bar over a symbol indicates the negation (or complement) of a variable. 

An underscore with a Boolean variable designates the positive logic convention; 
the absence of the underscore represents the more common negative logic convention. 

Logic Design 

. AND (conjunction or intersection) 

(The dot may be omitted.) 

+ (INCLUSIVE-)OR (disjunction or union) 

© EXCLUSIVE- OR 

implication 
Cl inclusion 

Boolean characterization 

(upper left line is input signal A; 
lower left line is input signal B; 
right line is output signal C) 

C = AB 
C = A + B 
C = AB 
C = A + B 
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gated pulse generator 


c = « a b 


C = Attg 

(No B input) 

S J 

flip-flop J where Sj, Rj, Tj are set, reset, and trigger inputs 
and J, J are assertion and negation outputs, respectively 


GENERAL THEORETICAL CONSIDERATIONS 

In generating a curve from any input point a to the next point b with slope /3 
at b (fig. 1), \j/ (the angle of b with respect to the presently generated position p) 

may be determined and eas- 
ily transformed to a new 
(primed) coordinate system 
(fig. 2) whose origin is at b 
and whose X-axis is alined 
with j3. The tangent angle <p' 
may then be determined from 
the simple relation 0' = 2t//' 
and easily transformed back 

to the original coordinates to F '9 ure L_ Interpolated curve through input data points a. b, c, d, e. 

guide the curve generation of 
one increment. The process 
may be iterated until point b 
is reached. (See appendix.) 

The relation <£' = 2xf/' 
is a special case of 

<t>' = Q^' (1) 
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Figure 2 Portion of figure 1 transformed to new 
coordinate system. 



which, for appropriate values of the parameter Q, is the equation of well-known geo- 
metric curves which go through the coordinate's origin at zero slope; for example, Q 
values of 1.5, 2.0, and 3.0 yield cardioids, circles, and lemniscates, respectively, where 
the size of these curves is determined by the initial point a. 

The subject controller's operation, which is based on the foregoing principles, is 
outlined and illustrated by an example in the next section. 


Computational Algorithms 

Point p of figure 1 is moved in finite steps from a to b along a path ending 
with slope /3 by use of the following algorithm: 

(1) Compute the direction angle \p of the vector pb. 

(2) Transform ip into expressed in a new coordinate system centered at b 
and alined with ft 

(3) Compute <p' = Qi p' . 

(4) Inversely transform <p' into <p, expressed in the original coordinate system. 

(5) Use <f> to predict the next increment of the curve, thereby establishing a 
new p. 

(6) Repeat steps (1) through (5) until p and b coincide. 

A step-by-step example of this simplified algorithm, in which it is assumed that 
(a) the coordinates of points a and b are, respectively, x a = 3, y a = 2 and 
x b = 10, y b = 12, (b) for x a S Xp < x b , /3 = 20° (for the computation of ft see p. 20), 
and (c) the desired curve is circular and hence Q = 2, is as follows: 


( 1 ) 


ip = tan“l 


y P - 7b\ 

Xp - X b j 


= tan 


-1 




55° (since initially point p 


coincides with point a) 


(2) \p' = \p - j3 (appendix) 
= 55° - 20° = 35° 


(3) <p' = Q\p' 

= 2(35°) = 70° 

(4) <p = + 0 

= 70° + 20° = 90° 

(5) Since $ = 90°, "increment" p one unit in the plus y-direction. 

(6) Since the new p does not yet coincide with b, return to step (1). 
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(1) * = tan-l(f^i§) = 52° 

* 

(2) xp' = 52° - 20° = 32° 

(3) <p' = 2(32°) = 64° 

(4) 0 = 64° + 20° = 84° 

(5) Since the motion of an incremental orthogonal two-axis system is restricted to 
some multiple of 45° and since 90° is the closest multiple to 84°, again increment p 
one unit in the plus y-direction. However, simply to ignore this 6° discrepancy and sim- 
ilar ones for the various increments of the path would result in the generation of an 
octagonal shape instead of the desired circular one. The problem is solved by integrating 
these discrepancies in accordance with the method developed in the section "Derivation 

of the Output Function” and summarized on page 17. 

(6) Return to step (1). 

The foregoing procedure is repeated until p and b coincide, at which time the 
entire process is repeated for the next segment* of the curve. In general, there will be 
a slope discontinuity at the beginning of each new segment of the curve; this condition is 
rectified by the method presented in the section "Merging Function Criteria." 

The algorithm is now refined with regard to the computation of xp. The relation 
xp = tan -1 (f(x,y)) is relatively difficult for machines to compute as seen in reference 4. 
Unfortunately, step (1) requires this computation for each increment of the curve. Alter- 
natively, however, the algorithm may be changed to utilize, instead, the more easily com- 
puted 6 xp, the finite change in xp. The refined algorithm is as follows: 

(1) Compute the direction angle xp of the vector pb. 

(2) Transform xp into xp'. 

(3) Compute <p' = Q xp' . 

(4) Transform cp' into cp • 

(5) Increment p in accordance with <p. 

(6) Compute bxp due to the movement of p. 

(7) Increment cp by Q bxp. 

(8) Repeat steps (5) through (8) until p and b coincide. 


The interpolated curve is composed of segments connecting adjacent data points. 
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This algorithm is advantageous in that the slow computations of steps (1) through (4) 
are performed once, with only the fast operations of steps (5) through (8) iterated for 
each increment of the curve. 

Step (7) is based on the difference equation 

*k+i = + Q 6 ^k 

Since, by the fundamental difference definition 

^k+1 = ^k + 6< ^k 

it is necessary to show only that 

60 k = Q 6l f/y. 

By the invariance property of the linear transformation (appendix), 

6<£ k = 60 k 

By the fundamental definition, 

^k = ^k+l - K 

From equation (1), 

6<£ k = Q^ k+1 - Q^k 
Again, by the fundamental definition, 

60 k = Q 6i^ k 

Finally, by the invariance property, 

6$ k = Q 6^ k 

Derivation of Output Function 

Step (5) of the refined algorithm requires that p be incremented M in accordance 
with The coarse 45° angular resolution inherent in the output of an incremental, 

orthogonal, two- axis system poses a problem, in this instance, since such a system must 
generate 0, a function of considerably finer resolution. 

The problem is first approached in a precise manner. From this approach a com- 
plicated partial solution is obtained. Then applied is an approximation which yields a 
complete and greatly simplified solution without, for most applications, excessively 
degrading the results. 
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Figure 3 shows a magnified portion 
of the output grid on an incremental device, 
function, and so forth. Intersections of 
grid lines represent stable or valid output 
locations, quantities, and so forth. 

Curve s represents the desired, com- 
puted, interpolated path or function where 
<fi is the angle of the vector or direc- 
tional tangent^ to the curve. For pur- 
poses of this discussion, the grid spacing 
is considered small enough with respect 
to the curvature of the path d$/ds to assui 
ments shown in figure 3. Other parameters 



Figure 3.- Magnified portion of incremental output grid. 

ne (p to be constant over the two incre- 
in figure 3 are defined as follows: 


0 desired motion angle necessary for coincidence with s at its next grid 

crossing 


6 C actual motion angle 


<9 f 


difference angle defined implicitly by 


e = B c + 0f (2) 

where 0 C an< * ma y be construed as coarse and fine components of 6 and are sub- 
ject to the two following constraints: 

0 ^ \o c = 0(mod 45°)J < 360° 

and 


-22.5° ^ 6 f < 22.5° 


p,r 


predicted and corrected error distances, respectively, between grid corners 
and s-curve — grid- line intersection 


k-l,k,k+l 


subscripts denoting order of the motion increments 


2The tangent is not the usual scalar tangent since the direction or sense in which 
the interpolated curve is generated affects <p - that is, if the sense of generation is 
reversed, then 4> is rotated by 180°. 
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The solution proceeds as follows: P^ + i> the error distance which would result if 
the motion continued in the direction of 8 C is calculated. If |Pjj + i| =0.5 incre- 
ment, then 


and 


If [Pjj+il >0*5 increment, then 


^c,k ^c,k-l 


y k+l “ p k+l 


0 c,k~ 0 c,k-l ±45 ° 


and 


r k+l " p k+l ± 1 

With the aid of figure 3, the procedure may be formulated for -45° ^ 8 C ^ = 45° as 

p k + l = y k + tan< ^k- tan0 c,k-l (3) 

r k+l = p k+l “ 1 


0 c,k " 0 c,k-l + 1 


where 


l[] is defined as the greatest integer in quantity within the brackets, and tan 8 Q is 
necessary, when 8 Q = ±45°. 

The foregoing formulation will be used in the derivation of equation (7); first, how- 
ever, it is necessary to establish the relationship between p and 8 £. Figure 3 shows 
that p may be expressed also as 


p k + l = tan % = tan ( e c,k + \k) 


(4) 


p k+l 


tan ^c,k + tan g f,k 
1 - tan 8 C> k tan 8f ^ 


(5) 
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For 0 c k = 0, equation (4) reduces to p k+1 = tan k . Hence, a change in 0 Q k should 
take place when 


0 f,k “ tan_1 ( p k+l) | > (| tan " 1 (±°-5)| = 26.6°) 

For 0 c k = +45°, equation (5) reduces to 

1 + tan k cos 0 f k + sin 0 ^ k 
p k+l l - tan 0fk cos 0f^ - sin 0f k 


cos 2 0 f k + 2 cos k sin k + sin 2 0 f k 

cos \ k - sin \k 


1 + 2 cos 0 f k sin k 
1-2 sin 2 k 


1 + sin 20 f k 
cos 26j k 


p k+l = sec 20 f,k + tan 20 f,k ( 6 ) 

Although the inverse 0 ^ k (P k+ j) is more complicated, it is easily seen by substitu- 
tion that a change in 0 Q k should occur when 

|0 f k | > (45° - 26.6° = 18.4°) 

that is, when 0^ k < -18.4°. 

Finally, for 0 C k = -45°, equation (5) reduces to 

- 1 + tan 0 f k 
1 + tan k 

Following the foregoing procedure results in 

p k+1 = tan 28 fk - sec 2e, k 

and, hence, a change in 0 Q k should occur when 0 f k > +18.4°. 

Even considering only these limited cases, the P k+ j function is relatively com- 
plicated. To assist in the application of a reasonable approximation, the function is dis- 
played, with the aid of figure 4, in figure 5. It is noted that to implement p accurately 
requires the use of one complicated function if 0 C - 0° and another complicated function 
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One output increment 



if 6 C - ±45°. The approximation shown in figure 5, however, greatly simplifies the 
problem since 

(1) It is linear and thus is simply p = C0f where C is the proportionality 
constant. 

(2) The same function is applied for both conditions of 9 C . 

(3) Its linearity enables the application of the following superposition property for 

p(0): p(#i) + p(0 2) = + &2)’ Note further that C vanishes when 6f is measured 

in octants (1 octant = 45°) and p in increments - that is, p = Of. 

(4) It avoids the problem of large tangents and, therefore, can be used throughout 
the entire 360°. 

The penalty for using this approximation is exhibited in the digital computer simu- 
lation of reference 4, appendix C, where the attempted difficult task of generating a circle 
through two points and one slope resulted in four sides slightly flattened. For normal 
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Figure 5.- Relation between actual and approximate error angle 0f and error distance p. 

curve interpolation tasks, however, this slight tendency to flatten is insignificant. The 
implementation, therefore, is based on the use of this approximation. 

In effect, then, the approximation establishes an equivalence relation between 
angular dimensions at p k of figure 3 and linear dimensions at the p k+ j increment 
crossing, so that the tangent functions are eliminated and thereby the useful angular 
domains are extended throughout the entire 360°. Implicit in the relationship p = 0 f is 
the corrected form y = 9 f (where the constraints of equation (2) have been applied to 
the latter dfj. Hence equation (3) can be written as 

e f,k = e f,k-l + ^k ' 0 c,k-l ^ 

where the constraints of equation (2) are not yet applied to k . By updating the 
subscripts 

0 f,k+l = 0 f,k + \ " 0 c,k 
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where, as assumed, Therefore, 

60 f,k = e f,k+l " 0 f,k = ^k " 0 c,k 
Integrating finitely (ref. 7) gives 

k-1 

%,k = 0 f,o + ^ 60 f,j 

i=o 

where j is a dummy variable. By letting the initial condition 0.j q = 0, 

k-1 

0 f,k = X ^ 

3=0 

where e-j = 6 0 f . = p. - 9 C i, the error between the tangent angle of curve s and the 
J 1 »3 3 

actual motion angle. Hence equation (2) can be written 

k-1 

0 k = 0 e,k + Y, € j 
3=0 

When e j exceeds the constraints of equation (2), 1 or more octants are transferred 

from ej to 0 C . Equation (9) is a most important result since step (5) of the refined 

algorithm can now be stated explicitly: 

(1) Determine 9 j ^ for each point p^ along the generated path by finitely inte- 
grating or summing the* error e j from initial €q = 0 to e^-1’ continually applying 
the constraints associated with equation (2). 

(2) Use equation (2) and associated constraints to establish 9 C ^ for each 
point p^. 

(3) Generate the path in the direction of 9 Q 

Thus, e j as previously defined represents the error introduced at the jth step 
because the actual direction of travel 9 C j differs from the computed direction py 
Whenever 10* . I, the finite integral of this error integrated from an input data point pg 

’ 1 ° r. 

to the previous generated point p^_ j, exceeds 22— , it applies to 0 £ ^ a 45° correction 
which, in general, is an overcorrection, in which case e j changes sign and 0j ^ builds 
up in the opposite sense. Hence the actual path (controlled by 0 C ) continually seeks and, 
in general, follows the computed path (controlled by p) to within 0.5 increment. In 
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terms of the CalComp digital incremental recorder model 565 that was used in this 
research, this 0.5 increment represents an error of 0.005 inch, a distance approximately 
as small as the plotted curve thickness. 

The beauty of this formulation is that it is easily implemented. 

Merging Function Criterion 

In generating a curve from point a to point b, the refined algorithm of page 10 
gives no consideration to the angle at which the curve originally entered a. In general, 
disregard of this angle results in derivative discontinuities at the various input data 
points of a curve. 

To obviate this problem, a merging function is derived which effects a smooth 
derivative transition from the curve entering some point to the curve leaving this point 
by limiting dtp/ ds, the curvature of the path in the vicinity of the point. The particular 
merging function criterion which is implemented may be considered to be circular since 
it limits dcf>/ds to the dp/ds of a circle whose diameter is equal to the distance 
between the present position p and the next data point. For the circle shown in fig- 
ure 6 d<p / ds is constant over the path s; hence, 

deft _ 2tt radians _ 2 radians <■ 2 radians / 1Q \ 

ds irdq increments dq increment “ dp increment ' 

If dp/ds is limited to 2/dq for an interpolated path, the merging function gives, in 
general, (1) a smooth transition, and (2) complete merging before p coincides with b. 

Although other merging criteria could have 
been used, this particular one is attractive since 
it is a simple function of dp, the square of which 
is already computed by the experimental hardware 
for other purposes. 

The merging function is shown in figure 7 as 
the limiting dp/ds plotted against dp. Limiting 
dp/ds is given in radians per increment and ir 
radians per increment where lit radian = n radians. 

The implementation described in the section 
"Operations of the Experimental Model" is based 
on an approximation to this function. 
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Figure 6.- Circular criterion for merging 
function. 



Distance to next point, dp, increments 

Figure 7.- Merging function. 

Computation of /3 

The development of the theory is based on the hypothesis that the slope angle (3 
of figure 1 is given. An examination of this hypothesis is now made. 

In general, a unique accurate slope angle /3 does not exist since the input points 
do not specify the curve in the immediate neighborhood of point b -r that is, in digital 
terms, the input points do not specify the increments immediately preceding and following 
point b. Therefore, the slope assigned to point b is a function of the particular 


interpolation scheme used. On the assumptions that sufficient input points are given to 
faithfully represent the function and that point b lies in the middle 50-percent region 
of the distance between a and c ^or, more precisely, i < X < 3, where X is the ratio 
of the distance ab to the distance acV then the simple first difference definition of 


/3 = tan" 


y a - y c 
x a - x c 


is, in general, adequate as is seen by examining the geometry of various three-point 
patterns. 


Equivalence of Angles 

Physically xp, (b, fi, and so forth, can range only over 2 it radians. Hence, 
physical significance is extended to those values outside of range by letting the equiva- 
lence relation 

mean xp^ = + 2 7 tK for some integer K (i.e., xp^ s ^(mod 2ir)). 

Principal Range of Angles 

With the aid of the foregoing equivalence relation, xp' in any range may be con- 
verted to its principal range defined as 

-v = xp' < + 7 r 

The important relation <p' = Qp' tacitly assumes that xp' is expressed in its principal 
range. However, it is shown in the following theorem that this requirement is really not 
a necessary condition when Q is an integer; that is, for the relation = Q xp' equiv- 

alent angles 4>' are obtained for equivalent angles xp ' when Q is an integer. 

Theorem: Under the relations 

( f) 1 = Qxp 1 

<P 2 ~ Q x I / 2 

xp 2 = xp ^( mod m) where m is any real number, 
if Q is an integer then 0 2 s $^(mod m). 

Proof: (1) xp 2 = i^(mod m), by hypothesis. 

(2) xp 2 = xp^ + mK, by the equivalence definition. 

(3) 0 2 = Q^2’ by hypothesis. 
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(4) <fi 2 = + mK )> steps (2) and (3). 

(5) ^2 = ^1 + ky hypothesis. 

(6) Q is an integer, by hypothesis. 

(7) KQ is an integer, since the product of two integers is an integer. 

(8) Therefore - 0^(mod m), by steps (5) and (7) and the equivalence 

definition. 

The converse is not true. However, an important fact concerning the converse is 
stated in the following corollary to the theorem. 

Corollary: Under the relations of the theorem, the equivalence of and 

is guaranteed only if Q is an integer. 

This corollary is a consequence of step (7) of the theorem proof; that is, although 
specific cases can be found in which KQ is an integer for a noninteger value of Q, an 
integer product is assured only if both K and Q are integers. 

These considerations — angular equivalence, principal range, and integer values 
of Q — are extremely important in terms of hardware implementation. Their signifi- 
cance will be briefly mentioned so that they can be more fully appreciated when the actual 
design of the experimental model is presented in a subsequent section. 

The angle equivalence relation 

1^2 = ^j(mod 27r) 

implies that the hardware need not distinguish angles which differ by integral multiples 
of 277. Therefore the associated counters, registers, and so forth, are required to 
handle precisely the modulus 277 range and, hence, bit positions of weight ^277 need 
not be provided even though the angle may exceed 2?7. 

The theorem essentially states that the complication associated with the principal 
range constraint can be avoided if one is satisfied with interpolated curves involving only 
integer values of Q - that is, Q = 1, the linear case, Q = 2, the circular case, and so 
forth. Utilizing only integer values of Q would greatly simplify the system control 
problem and could eliminate the binary rate multiplier in the design of the experimental 
model. The importance of noninteger Q in a practical operational device can best be 
established through operational experience. The more complicated noninteger Q capa- 
bility is therefore retained in the model design so that some measure of its importance 
may be acquired through use of the experimental device. 
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COMPUTATION OF ip 


x 


Perhaps the most important single 
function to be performed in the interpo- 
lation scheme is the computation of 

* = tan' 1 p- 
Ax 

or, if one of the points is at the origin 
as in figure 8, 

ip = tan" * — 
x 

Thus, with a suitable arc tangent algo- 
rithm, ip can be computed absolutely. 

Such a computation with digital devices, 
however, is complicated and time con- 
suming relative to the more common 

algebraic manipulations. (See ref. 4.) Fortunately, in the present application an absolute 
computation is not, in general, required or desired because, once \p has been estab- 
lished for two data points, a and b, it is necessary only to add the incremental change 
in ip as the curve "steps" (with perhaps 100 increments) from point a to point b. 

To determine the incremental change dip in \p( dx,dy), proceed as follows: 



Figure 8.- Effect of fix and 6y on <i. 


tan \p = | 

d(tan ip) = d(|) 

o , , , x dy - y dx 
sec *ip dip = ^ 

x“ 


X 2 + y 2 X dy - y dx 

— OXy - 


and, therefore, 


dip = 


x dy - y dx 
x 2 + y 2 


(ID 


Thus, for each infinitesimal change in x and/or y, the resulting change dip can be 
computed. Furthermore, this change can be integrated over some path to obtain the total 
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angle ip . Independence of path, singularities, and other considerations of rigor are 
shown in advanced calculus texts (e.g., ref. 8). 

Unfortunately, the present application involves finite changes, to be denoted by 6x, 
6y, &\]s, and so forth, instead of infinitesimal changes. To combat the truncation errors 

caused by the finite integration process, a finite differential or difference is derived as 
follows: 

| = tan ^ 

and, hence, 


— + — = tan(i// + 5i//) = 
x + 6x 


tan ^ + tan 6ij/ 
1 - tan i// tan 6 \f/ 


| + tan 6)// _ y + x tan 
1 - | tan 6\p ~ x “ y tan 6^ 


Cross multiplying gives 

xy + x 6y - (y + 6y)y tan = xy + y 6x + (x + 6x)x tan 6i// 


and, hence, 


tan 6i//(x(x + 6x) + y(y + 6y)) = x 6y - y 6x 


When is small, 


6 ip ~ tan 6 ^ = 


x 6y - y 6x 
x(x + 6x) + y(y + 6y) 


( 12 ) 


By applying lim tan = 6i(/ and lim (x + 6x) = x, equation (12) reduces to 
6x— 0 5x— 0 

6y— 0 

x 6y - y 6x 

6i// = - — which, as might be expected, is the infinitesimal case. 

x 2 + y 2 

Equation (12) can also be written as 


x k x k + i + y^k + i 


(13) 


which shows that in the denominator of the finite case, both the present (kth) value and 
the next ((k + l)th) value of the independent variables are used. 

Equation (13) forms the basis of the experimental model design. Its implementa- 
tion is best understood by first considering equation (11) in terms of some of the better 
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known digital differential analyzer (DDA) techniques (refs. 9 and 10). Accordingly, the 
notation of figure 9 is interpreted as follows: 

dy pulse representing a small, constant increment of y(x), which is counted into 

(integrated by) register Y 

dx pulse representing a small, constant increment of x, which causes y(x) to 

be accumulated into an internal register R (fig. 10) 

dz pulse representing a small, constant increment of function z (dz is 

emitted when R overflows) 

Function z may therefore be obtained by counting the dz pulses - that is, 

pb pb 

z = Ja dz = Ja ydx 

This rectangular integration process is displayed in figure 11. 

One DDA approach to the implementation of equation (11) is the servo technique 
which follows. Rearranging this equation gives 

0=xdy-ydx- (x 2 + y 2 )d \j/ (14) 

In general, if there is a change in x and/ or 
y, equation (14) is no longer equal to zero but 
is equal to some unbalanced error a. This 
error can be nulled (i.e., caused to approach 

zero) by the generation of proper polarity Figure 9.- Conventional DDA integrator notation, 

dt// increments. Equation (14) can therefore 
be written 

dt^ (15) 

If the basic integrator of figure 9 and a few 
additional symbols are used, the configura- 
tion of figure 12 is evolved as the imple- 
mentation of equation (15). The additional 
notation shown in figure 12 is as follows: 

dy 

Figure 10.- Mechanization of figure 9. 


a = x dy 


J 2y dy 
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Figure 12.- Symbolic implementation of equation (15). 

a bidirectional counter called a summer whose input is time shared by three 
sources which appear as terms in equation ( 15 ) 

indicates that y dx and (x^ + y^)di/y are subtracted; its function may be 
considered that of reversing the counting direction of o 




pulse generator controlled by the output of a (The tri valued output (+, -, 0) 
causes the generation of +\p increments, increments, and no 
increments, respectively.) 

x2 signifies that the input is multiplied by 2 as it enters the counter (If the 

counter is binary, then this multiplication is accomplished simply by 
connecting the input to the next most significant stage.) 

The additional input to the interior integrator is time shared with its normal input. 

With further extensions of the notation, equation (13) may be represented by fig- 
ure 13. The only new symbol which appears is the double-headed integrator which emits, 
in addition to the usual kth output at one head, both kth and (k + l)th values at the second 
head. Note from the figure that the top input to the integrator goes with the top output 
and the bottom input goes with the bottom output. 

That the second head can satisfy equation (13) is seen as follows. At the end of the 
kth computation, assume that 

D = X 2 + yj (16) 

Assume that x changes 1 increment - that is, 6x = ±1. Let this change cause, at 
first, only the kth value to be emitted from the second head. Then, 

D = x k + y k + x k 6x 

- M** + Sx ) + y k 

■Vki-i + rk 
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Since there was no change in y, y k+ j = Y k and hence y k = y k Y k+1 - Therefore, 

D - ’VW + Vk*i 

as required by equation (13). At the end of the computation, the remaining (k + l)th 
value is emitted so that 

D = x k x k + i + yk y k + i + x k + i 6x 

Again, since there was no change in y, y k = y k+1 and hence y k y k+ j = y k+1 - There- 
fore, at the completion of the (k + l)th computation, 

D = x k + i( x k + 6x ) + yLi 
= x k+i x k + i + yLi 
= x Li + yk + i 

a result compatible with the assumption of equation (16). An analogous situation occurs, 
of course, for a change in y. 

Figure 11 shows that the quantity represented by one dz increment is relatively 
large. This large resolution of dz is necessary to insure that no more than one dz 

pulse is emitted for any integration step. Since the quantity represented by dz is 
C y max dx 

y dx, \ dz can be in error by ± . To eliminate this error and to facili- 

J max 7 J 2 

tate the design of the double-headed integrator, the following significant breaks with 
convention are incorporated in the implementation of figure 13: 

(1) The integrators contain no R accumulators. 

(2) The integrators contain input accumulators instead of input counters; this 
change permits, of course, the acceptance of quantized data instead of the conventional 
single pulses. 

(3) The summer labeled a in figure 12 accepts quantized information instead of 
the conventional single pulses. 

Thus, in the conventional diagram of figure 12, all interconnecting lines "carry” 
single pulses, whereas in figure 13 all interconnecting lines, except those labeled 6x, 

6y, and "carry" quantized information. 

In order to check the errors associated with equation (13) as i{/ is integrated 
over various paths, a digital computer program was written. Results obtained from use 
of this program indicated that, in general, the accuracy was well within the requirements 
of the present application. For program details, see reference 4. 
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OPERATIONS OF THE EXPERIMENTAL MODEL 


General Description of Model 

As indicated in the system block diagram of figure 14 and the photograph of 
figure 15, the input to the experimental model is paper tape, coded in the format of 
table I by a Flexowriter. The model's output device is a standard CalComp digital 
incremental recorder. 



Figure 14.- Simplified diagram of the experimental model. 
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Figure 15. 4 * Front view of the experimental model. L-65-2370 




TABLE I.- TAPE FORMAT 


—-Tape motion 

| SXXXXASXXXXRSXxj j XXRE 


x-data y-data 



Nota- 

Printed 

Coded as 

Detection 

Character 

tion 

symbol 

X O C 8 4 2 1 

code 

End of data 

E 

E 

* * * * * 

04 

Carriage return 

R 

None 

* * * * * 

84 

Space 

A 

None 

* 

OC124 



i 

< 

+ 

* * * 

oc 

Sign digit 

S 

( 


* 

xc 




T) 

* 





1 

* 





2 

* 





3 

* * * 





4 

* 


BCD digit 

X 


5 

* * * 





6 

* * * 





7 

* * * 





8 

* 





.9 

* * * 



The controller is referred to as INTROL, an abbreviation for higher order interpo- 
lation controller and, for the model, is totally contained within the relay rack in figure 15. 
It is implemented with the Computer Control Company 1-Mc S-PAC digital modules in 
which up to 28 logic cards are inserted in each of the 5^- inch S-BLOCs seen in the 
photograph. 

In general, the mission of INTROL is to accept input data points (similar to those 
of fig. 1) and direct the motion of the recorder along a smooth path (similar to the path in 
the fig.) through the points. The hardware which performs this operation is diagramed 
in figure 16. 

The box labeled "input data" at the top of the figure represents the tape reader or 
manual data entry switches which read input data into RC. When a new data point is 
required, as indicated when D = 0, the point is gated by block Gj into RC and the orig- 
inal contents of RC is shifted into RB and those of RB into RA. Switch 1 signifies 
that A1 can read the contents of any of the three registers as required for computing 
first differences necessary in the computation of arc tangents. 
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Figure 16.- Functional diagram of experimental model. 

The block labeled "i^-computer" represents the computer previously shown in fig- 
ure 13 with the addition of a binary rate multiplier; it should be viewed now simply as a 
block in which 6x and 6y increments are applied and the resulting angular change A<£ 
is emitted. The D output, where D is the square of the distance to the next point, is 
used not only in the computation of A <p but also by the merging logic to establish the 
instantaneous merging rate. Blocks C2, C3, and C4 are counters; C4 counts in accord- 
ance with the constraints associated with equation (2). The direction of 9 C , whose reso- 
lution is 45°, is interpreted by the theta decoder block in terms of 6x and 6y compo- 
nents. These 6x and 6y increments drive the CalComp recorder and are fed back to 
the ^-computer for determining their effect on <fi and D. Block G2 is used to preset 
(C3) to -6 C as an initialization, where (C3) means the quantity in C3. 
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Simplified Operations of Model 

An operations algorithm of the experimental model is displayed in figure 17 in the 
form of a flow chart. Although considerable detail is omitted, the chart expands the 
refined algorithm given in the section "General Theoretical Considerations" to include 
the operations involved with the output and merging functions subsequently developed in 
that section. 



Figure 17.- Simplified operations flow chart of experimental model. 
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The operations algorithm for the hardware of figure 16, which is indexed with the 
numbered boxes of the flow chart of figure 17, is as follows (assume that the curve in 
fig. 1 has been generated up to point a and that signal D therefore has been reduced 
to zero): 


(T) Shift data from RB to RA and from RC to RB and read a new point of 
the curve into RC. (The points are plane Cartesian and hence the registers and input 
accumulator A1 are actually implemented in x-y pairs.) 


(D Reset Cl and C2 and A1 and A3; preset (C3) to -6c- 

(3) Switch A1 first to RA and then to RC to accumulate (Alx) = x a - x c 
and (Aly) = y a - y c - 


(4) Integrate, with Cl which is internal to the ^-computer, 5\p over some path 
to obtain -/3 = -tan - ^ Simultaneously, multiply the same 6i// pulses by 

- ^ in the binary rate multiplier and then integrate with C3. (Thus, 

(C3) = -j3(Q - 1) - 6 C expressed in n radians. ^ That the binary rate multiplier gives 
the desired result can be seen in step (6).) 


(5) After resetting Al, accumulate (Alx) = x a - x b and (Aly) = y a - y^. 

(§) Determine, with the ^/-computer, tan - ^ by integrating, without reset, 


the associated 5xj/ pulses in Cl, making (Cl) = -/3 + ^ - 4 /' . Simultaneously, multiply 
these dijy pulses by Q/tt in the binary rate multiplier and integrate the product in C3, 
making 


(C3) = Ofr - j3(Q - 1) - 6 C 

- Qi'P - 0) + fi - 9 C 

- Qip' + ~ &c 

= 0 T + / 3 - 6q 


- (f) - 6 C 

= ^ - e c0 

where 0 C 0 = initially (i.e., the slope associated with point a) and 6 C( p ~ <f> after 
merging is complete, at which time (C3) = 0. (Cl, which contains ij /' , as is shown in a 
subsequent flow chart, is used to constrain \j/' to its principal range (p. 20).) 

^ The conversion from radians to n radians simplifies the theta decode logic. It 
also keeps the multiplying factor Q/u less than unity (for Q < tt) as required by the 
binary rate multiplier. 
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(7) For merging, reduce |C3 J by A0 C( ^ = f ((A3), S[C3]), where (A3) = D is, 
again, the measure of the distance to the next point. Simultaneously, change (C2) by the 
amount of A 9 C( p. (Note that (C2) = 6 C( p - 9 C = e, where, as previously explained, 9 C( p 
merges in amounts of A 9 C( p from 0 C to cf> and hence e changes from zero to 
<p - 9 C , the difference between the computed direction angle 4> and the actual 45° reso- 
lution angle of motion 0 c .j 

(D With C4, integrate the error e contained in C2 as required by equation (9). 
If a change in the coarse portion of C4, 6(C4<j) = 60c> occurs, feed this change back 
into C2 so as to update its Q c . 

(9) Since, during the first traversion of the major loop, (A3) 4- 0, proceed. 

(10) With the theta decoder, resolve 9 C into the two orthogonal increments 6x 
and 6y . 

(y) Using the bipolar pulses representing the 5x and 6y increments, step the 
CalComp recorder and increment the ^-computer (so that it updates xp and D). 

(12) Since merging would not normally be complete (i.e., (C3) =£ 0), proceed. 

(l^) Change (C3) by A <p = Q Axp, where Axp = Z5xp for a given 6x or 6y 
increment, and then enter the major loop at step (7). 

(l|) If merging is complete in step @) (i.e., (C3) = 0, change (C2) by QAxp 

and then enter the major loop at step (8). 

Eventually the test at step (9) will show (A3) = 0 at which time the entire operation from 
step (T) is repeated. 


Detailed Operations of Model 

The complete operation of the experimental model consists of the following 
programs: 

(1) ^/-computer subroutine 

(2) main program 

(3) input program 

The ^/-computer subroutine is called only by the main program. The input program 
runs simultaneously with the main program except for an interlock flip-flop designated G. 
The main program requests a new data point from the input program by setting G. After 
making the new point available, the input program resets G. After using the point, the 
main program again sets G, and so forth. 
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^-computer subroutine .- The i//-computer, which now includes a binary rate multi- 
plier, is treated by the main program as a subroutine (fig. 18). It is called with the 
statement: 

Call xj /- comp (6x,6y); A <p) where 6x,6y increments are input arguments and 
A<£ = 26 cp, an output argument, is the total change in <p resulting from a 6x and/or 
6y increment. 

The quantity A 0, which appears also in figure 16, is a convenience and does not 
actually exist — that is, there is no register which contains A<)f) = 26 <p. Instead, C2 
and C3 actually utilize the individual 6 <p pulses. The actual representation, however, 
would require a transfer between the main program and the subroutine for every 6 <p 
pulse iteration. Fortunately, this complication is eliminated by the fictitious A <f> 
accumulator. 

Although it is not explicit in the subroutine call statement, the main and subroutine 
programs can each sample the other’s quantities, registers, and so forth. 





The i//- computer subroutine flow chart of figure 19 should be used in conjunction 
with the i//- computer block diagram of figure 18 which is almost the same as that of fig- 
ure 13. In the following explanation of the subroutine flow chart, the numbered explana- 
tion steps are indexed with the numerals at the flow chart blocks: 
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(T) Reset the fictitious A<p accumulator. 

(2) If J 6y | = 1 and |6x| = 0, proceed. 

(3) Reset flip-flop V. (Flip-flop V then causes computer to operate on 6y.) 

(?) Accumulate in A3 the product of (A2y) and 5y. 

Accumulate in A4 the product of (A2x) and 6y. 

Accumulate in A2y the 5y increment. 

(§) Reset 6y. 

(§) If, in step (2), |6y | = 0 and |6x| = 1, proceed with steps (6), (7), and (§) in an 

analogous manner to steps (3), (4), and (5). 

(D If, in step |6y| = 1 and j 6x J = 1, proceed with step (§). If V is reset, 
go to step (4); if V is set, go to step (7). 

@ If (A4) =* 0, proceed. 

(fi) Assign the sign of (A4) to 

(f2) Accumulate in A4 the negative of the product of and (A3). 

(R) If flip-flop 0 is set, proceed. (0 = 1 signifies that angle 0 is being 

computed.) 

(£?) "Decrement" (Cl) by 6 if/ and set F = -(Q - 1 )/t r. (See step (5) on p. 32.) 

(f^) If, in step (£3), 0 = 0, which signifies that \p is being computed, increment 

(Cl) by 6 \p and set F = Q/ir. (See step (6) on p. 32.) 

(f(^ Accumulate in the fictitious A (p accumulator the product F Go to 
step (fiJ). 

(T?) If (A4) = 0, proceed with step (f?). If V = 1, proceed. 

(f§) Compute AA3 = x^ + j6x. 

(I§) Accumulate AA3 into (A3). Go to step (2). 

(gft) If, in step (f?), V = 0, compute AA3 = y k+1 5y. Go to step (T§) . 

(0) If |6y| = 0 and 1 6x | = 0, return to the main program. 
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Figure 19.- (//-computer subroutine flow chart. 
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Main program.- An appreciation of the procedure used in integrating is helpful 
in understanding the main program flow chart. Recall from equation (12) that it is 
important to keep the change in i// small for incremental changes in x and y. Also, 
the angle in figure 20 may be integrated by jumping from the origin to point m and 
then integrating to point a as was done in reference 4. This jumping avoids the singu- 
larity (ref. 8) at the origin, saves integration time, and is permissible since the path 
from 0 to m has no effect on i//. In terms of the hardware of the experimental model, 
this jumping presents a problem since (A3) of the t//- computer must be preset to the 
product xj[, a quantity which, in general, is unknown. Hence, the design of the experi- 
mental model is based on the following integration procedure: 

(1) Jump exactly 12 increments in the negative x-direction. 

(2) Integrate 12 increments in y in the direction of the data point (i.e., point a 
in fig. 20). 

(3) Integrate along the x-direction toward the point until its x-coordinate is reached. 

(4) Integrate along the y-direction toward the point until its y-coordinate (and hence 
the point) is reached. 

Such a procedure avoids the singularity at the origin, permits presetting (A3) 
always to 12^ = 144, and keeps the changes in xj/ small. 
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Figure 20.- Integration of angle ip. 
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To detect when the coordinates of the point are reached as required for steps © 
and (3) in the integration procedure, (Alx) is decremented for each step in the 
x-direction and (Aly) is decremented for each step in the y-direction. A coordinate 
of the point is reached therefore when the associated accumulator reaches zero. 

The main program flow chart of figure 21 should be used in conjunction with the 
experimental model hardware of figure 16. In the following explanation of the main pro- 
gram flow chart, the numbered explanation steps are indexed with the numerals appearing 
at the flow chart blocks: 

(1) Select the desired type of trajectory with switch Q. Set recorder to initial 
position, lower pen, insert paper tape, and so forth; push start. 

(2) Reset C3 and flip-flop L; preset U to 1 and set flip-flop G. (Flip-flop G 
interlocks the main and input programs which otherwise operate independently.) 

©If G = 1, proceed. 

(?) If end of data character E has not been read on tape, go to step (3). 

© If, in step ( 5 ), G = 0, which signifies RA, RB, and RC are filled, proceed 
with step ©. If U < 3, proceed. 

© Increment U. 

© Set G; go to step ©. 

© If, in step ©, U = 3, reset Cl and A1 and set flip-flop (3. 

© Preset (Alx) to 12. (See step © of integration procedure.) Accumulate 
(RAx) and -(RCx) so that (Alx) = x a - x c + 12. 

( 1 ( 5 ) Accumulate (RAy) and -(RCy) so that (Aly) = y a - yc- 

© Reset A2, A3, and A4. 

(O) Preset (A2x) to -12 and (A3) to x^ = 12^ = 144. (A2x and A3 are a 
part of the ^-computer as shown in fig. 19. For an explanation of presetting to -12, 
see step © of integration precedure.) 

(O) If (Aly) = ya “ Yc is negative, proceed. 

(I?) Set 6y = -1 and 6x = 0. Go to step 

(£5) Set J to 1. Go to step @. 

(Ijj) If, in step (f§), (Aly) is positive, set 6y = +1 and 

@ Decrement (Aly) by 6y and call \f /~ comp (6x,6y; 
begins the integration of (3.) 


6x = 0. Go to step (Ij}) . 
A0). (The i//- computer 
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© Increment (C3) by the ^/-computed A$ = SS<£. 

0 Repeat steps (17) and 0 eleven more times. (See step (2) of integration 
procedure.) 

(0 If (Alx) is negative, proceed. 

(§1) Set 6x = -1 and 6y = 0. 

(0 Decrement (Alx) by 6x and call i//-comp (5x,6y; A<p). Go to step 0). 

0 If, in step (0) , (Alx) is zero, set 6x = +1 and 6y = 0. Go to step (0 . 

(0 Increment (C3) by A 4>. 

(0 through @ Treat in analogous manner to steps (0) through (0 . (Integra- 
tion of /3 is complete.) 

(0) If, in step (0) , (Aly) is zero, proceed with step (0) . If 0=1, proceed. 

(0 Reset flip-flop 0, which signifies that xp is to be computed. Although 
redundant, reset Al. 

(0) Preset (Alx) to -12. Accumulate (RAx) and -(RBx) so that 
(Alx) = x a - Xfo + 12. 

(0) Accumulate (RAy) and -(RBy) so that (Aly) = y a - yjj. Go to step (O). 

(0) If, in step 0), flip-flop 0 is reset, set flip-flop G. (Cl now contains 
_j 3 + xp = xp\ the transformed angle.) 

0 If | Cl | > 7T, which signifies that xp 1 is not in its principal range, proceed. 

0 Set J to 402. (Note that since the weight of one 6 xp pulse is 2"® radian, 

402 6 xp pulses represent 2i r radians. Therefore, 402 pulses (or integral multiples 
thereof) are used to bring xp ' within its principal range. <p = (C3) is changed 
accordingly.) 

(0 If (Cl) is negative, proceed. 

0 Give 5xp a negative sign. 

(0 If J > 0, go to step 0 ; otherwise, go to step (0 . 

0 If, in step 0, (Cl) is positive, give 6xp a positive sign. Goto step 0. 

(0 If, in step (0, J > 0, add 6 xp to (Cl) and F 6xp to (C3). 

0) Decrement J. Go to step (0 . 

0 If, in step 0, |Cl| ^ ^ r, proceed with step 0. If U 13, proceed. 

(xp' = (Cl) has been computed and confined to its principal range and <f> is contained 
in C3.) 
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0 ) Set the coarse portion of (C4) equal to the coarse portion of (C3). Reset 
(C3 c ), leaving (C3) = <j) - 6 C . (Before the first portion of the curve is plotted, a reason- 
able initialization should be assigned to (C4 c ), the slope of the curve at the first point.) 

0 ) Reset C2. 

0 >) If (A3) > 0, proceed. 

0 ) Decode (C4 c ) into its 6x and 6y components. 

(0) Increment the recorder with 6x and 6y. 

(0) Increment the ^/-computer with 6x and 6y. 

0 If (C3) * 0, proceed. 

0 Increment (C3) by A<p = Z6<p from the ^-computer. 

( 0 ) Increment (C2) by A 9 C $ = f((A3),S[C3]), the merging quantity obtained from 
the ^/-computer. Decrement (C3) by A0 C( ^. 

( 0 ) If (C3) changes sign, proceed. 

( 0 ) If (C3) # 0, proceed. 

(0) Decrement (C2) and increment (C3) by tt2~^ with a sign identical to 
S [C3] . Go to step ( 0 ) . 

( 0 ) If, in step ( 0 ) , (C3) = 0, proceed with step ( 0 ) . 

If, in step ( 0 ), (C3) changes sign, proceed with step 0 ). 

Increment (C4) by (C2) and decrement (C2) by the amount of the change 
in (C4 c ). Go to step 0 ). 

( 0 ) If, in step ( 0 ), (C3) = 0, increment (C2) by A 0. Go to step 0 ). 

0 ) If, in step 0 ), (A3) = 0, preset (C3) = -(C4 c ). Go to step ( 5 ) to interpolate 
the next segment of the curve. 

0 ) If, in step (?), end of data character has been read on tape, proceed with 
step 0 . If Z = 0, proceed. 

0 Set Z = 1. Go to step ( 8 ). 

0 If, in step 0 , Z = 1, stop the operation since the complete curve is now 
drawn. (Note that step 0 ) is equivalent to reading the last point twice.) 
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Figure 21.- Continued. 
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Figure 2L- Continued. 
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Input program.- The input program as previously indicated operates independently 
of the main program except for the interlock flip-flop G which, when set by the main 
program, starts the input program. The flow chart of figure 22 should be used in con- 
junction with the input diagram of figure 23 and the tape format of table I. The following 
flow-chart explanation steps are indexed with the numerals appearing in the blocks on the 
flow chart: 

(T) If G = 0, wait; if G = 1, which signifies that the main program has requested 
data, proceed. 

(2) If model is operating in manual entry mode (i.e., not in the paper tape input 
mode), proceed. 



Figure 22.- Input program flow chart of model. 
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(3) Set data point into manual switches. 

(3) Set J = 0. 

(5) Shift RA, RB, and RC in linear mode. 

(S) If J < 12, proceed. 

(7) Increment J. 

(8) If, in step (6), J = 12, transfer in parallel the data entry switch settings into 
RC (both x and y portions). 

d) If, in step (§), the model is operating in the tape mode, reset tape input register 
T and start tape reader. 

@ Read paper tape character. 

(Ti) If first tape character is a negative sign, proceed. 

(O) Set flip-flop S. 

@ Set L. Go to step @ . 

(O) If tape character of step (O) is a positive sign, reset S. Go to step (O) . 

@ Set J = 0. 

(fB) Shift (one step) tape input register in the circular mode. 

(T^) Repeat step (fB) eleven more times. 

(TB) Store BCD digit in least significant end of tape input register RT. Go to 
step (IB) . (Repeat steps @1 through (IB) for next three characters which are also BCD 
digits. This procedure reverses the digits for the ensuing BIDEC operation.) 

(IB) If, in step (O) , character is A, set flip-flop X which connects RT to RCx. 

(£5) Reset L which changes RT to linear shift mode. 

@ Set J = 0. 

(SB) If quantity is negative (i.e., S = 1), proceed. 

(SB) Route quantity through two's complementing circuit. 

(0 If, in step (SB), quantity is positive (i.e., S = 0), shift (one step) RT, RA, 

RB, and RC. 

@ For any decade D containing a number greater than 7, proceed. 

(SB) Decrement (D) by 3. 

(£7 ) In step (SB), for any decade containing a number less than 8, proceed with 
step (SB) . Repeat steps (0) through (BB) eleven more times. 
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(0) If X = 1, go to step (R5) . 

(S§) If tape character of step (0) is a carriage return, reset X. Go to step (£d) . 
(0) If, in step (^|), X = 0, set G = 0, which indicates to main program that new 
data have been entered. Go to step (I). 

(§1) If, in step (O) , character is E, which signifies end of tape, stop. 

The main program, xf /- computer subroutine, and input program form the basis for 
the design of the experimental model. 



Figure 23.- Data entry diagram. 

DESIGN OF THE EXPERIMENTAL MODEL 

The discussion of the design which follows points out how the subroutine, main, and 
input programs are implemented in the experimental model and proceeds in this same 
order. The logical design, in general, is presented directly in terms of NAND modules 
without the AND/OR design which sometimes precedes a NAND/NOR design. No attempt 
is made to show pin connections, duplications of circuitry, complete integration and con- 
trol circuits, and so forth; instead, included is only that portion of the design philosophy 
which is necessary for easy extrapolation to the complete wiring design. 

Figures 24 and 25 show the logic portion of the experimental model. By use of the 
rack layout of figure 26 the logic for the various functions performed by INTROL can be 
identified. In figure 24 is a temporary panel which corresponds to location S-BLOC 4 and 
which will be replaced by an operation flow indicator panel when INTROL becomes fully 
automatic. 
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Design of the I/'- Computer 


The ^/-computer occupies locations S-BLOCs 6 and 7, with its binary rate multi- 
plier in S-BLOC 8. Accumulators comprise a major portion of the computer. The 
accumulators are similar, but they do have significant differences. Therefore, a typical 
accumulator is described after which the peculiarities of each are pointed out. 

Typical accumulator .- In general, parallel adders are considered to be faster than 
serial adders. In some instances, however, the well-designed parallel adder is actually 
slower than the serial one that is implemented with the same speed logic. For example, 
the serial add time per bit for 1 megacycle logic is easily 1 /isec, whereas in a simple 
parallel adder 2.2 psec was required per stage to provide for carries. Obviously 
unless some provision is made perhaps to eliminate this time per stage for those stages 
which do not actually carry, the parallel adder is slower; to make this provision, however, 
greatly complicates (ref. 6) the circuit. Therefore, INTROL utilizes the serial natural 
binary type of accumulators. 


Figure 27 shows a typical serial 
accumulator in which the quantity B 
is added (or subtracted if N = 1) to 
(or from) quantity A and the result, 
replacing the old value, is stored in 
register A. The carry (or borrow) 
output P^ + j is delayed one time 
cycle r and fed back to the input. 

More typical for INTROL is the 
accumulator of figure 28 in which the 
B register is not an internal part of 
the accumulator and the carry (or 
borrow) delay is effected by the prop- 
agation flip-flop P. 


Table II is a truth table for the 
typical accumulator. Note that the sum 
A^ + j for addition where N = 0 = + is 
the same as for the difference A k+ j 
for N = 1 = -. Hence, for either addi- 
tion or subtraction, 


A k+1 = ABP + ABP + ABP + ABP (17) 
A k+1 = A ( BP + BP ) + A ( BP + BP ) 



Figure 27.- Typical serial accumulator. 



Figure 28.- Typical INTROL accumulator. 
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TABLE n.- TRUTH TABLE FOR TYPICAL ACCUMULATOR 


Input 

Output 

N 

A k 

B k 

p k 

p k +l 

A k +1 

0 = + 

0 

0 

0 

0 

0 

0 

0 

0 

1 — / 

— «- 0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

1 

1 

0 

0 

1 

0 

0 

0 

1 

0 

1 

0 

1 

1 

0 

0 

1 

1 

0 — < 

2— 1 

0 

0 

1 

1 

1 

1 

1 

1 H * 

0 

0 

0 

0 

0 

1 

0 

0 

1 

1 

1 

1 

0 

1 

0 — < 

t — «- 1 

1 

1 

0 

1 

1 

1 

0 

1 

1 

0 

0 

0 

1 

1 

1 

0 

1 — / 

3—0 

0 

1 

1 

1 

0 

0 

0 

1 

1 

1 

1 

1 

1 


where the kth subscript is implicit. Let 

H = BP + BP = B © P (18) 

and then 

A k+1 = AH + AH =A©H=A©B©P (19) 

Since EXCLUSIVE-OR logic was not readily available, 
A k+ j was implemented in the form of equation (17). 

The truth table shows that the propagation signal 
undergoes only 4 transitions from P k to P k+ j out 
of the possible 16. This fact leads to the use of flip- 
flop P and its simple logic which was obtained from 
figure 29 as 

Sp = NAB + NAB = B(NA + NA) = b(n © a) = BW (20) 
and 

Rp = NAB + NAB = b(nA + NA) = b(n © a) = BW (21) 
where W = N © A. 

The NAND implementation of these equations is 
shown in figure 30. 



Figure 29.- Map of propagation function. 
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Figure 30.- Implementation of propagation function. 


Accumulator 3 .- The accumulator of figure 28 closely represents A3 in which shift 
register A is conventional (ref. 4) and comprises 24 stages. Input B is time- shared 
ORed with A2x and A2y outputs. 

Accumulato r 2.- Accumulator 2, in which the x- and y-portions are identical, is 
required to perform only as a bidirectional counter in the sense that its input consists 
only of single plus or minus pulses. However, since A2 must serve also as the B 
register for A3 and A4, it must be capable of shifting. Also, although not required in 
the present design, both A^ and A k±l outputs are available simultaneously. There- 
fore, A2 is implemented as an accumulator. Its design differs from that of the typical 
INTROL accumulator only in that the B^ input (and associated logic) is eliminated. The 
counting is accomplished by setting the propagation flip-flop P and then by accumulating. 
Therefore, for A2, where B = 0, input equations (17), (20), and (21), respectively, 
reduce to 

A W = AP + AP 

Sp = J 6x | (for A2x) 

and 

Rp = NA + NA 

It is seen in figure 30 that, since B^ = 0, only the reset portion of the propagation func- 
tion is required for A2. 

The twelve stages provided in register A enable a capacity of sign and 2* 1 counts 
which are equivalent to a maximum distance of ±20.48 inches on the CalComp recorder. 
Since (A3) = x 2 + y 2 , the maximum value of 

(A3) =(2 11 ) 2 + (2 11 ) 2 
= 2 22 + 2 22 
= 2 23 

Therefore, 23 natural binary stages plus a sign stage are provided. (The sign bit in this 
application actually serves only as a check since it should never indicate negative.) 
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Accumulator 4.- Accumulator 4 is a typical INTROL accumulator except that the 
input is time shared by three sources. For the nulling operation (eq. (14)), the nulling 

direction is opposite to the sign of (A4) (i.e., N(A4] = S [A4] ) and shall be considered 
complete when (A4) changes sign. This change is detected by observing the propaga- 
tion flip-flop P after each accumulation; a change is implied only when P is found in 
the set state. 

In order to establish the sign and bit weight of A4, the following symbol definitions 
are used: 

Lpj length of (i.e., number of bit positions or stages in) RJ where J is any 

alphabetical character 

W5 U weight associated with least change in u where u is any variable or (AJ) 

To accommodate A3, ^ 24. Since A4 is nulled after each input cycle, it 

need be no larger than A3. It is recalled that Wg^3 = 2®. A resolution requirement 
of <1° was imposed on the experimental model; hence, was made 2“® » 0.9°. 

Therefore, to accommodate A3, Wg^ = Wg^gWg^ = 2® • 2"® = 2“®. To accommodate 

the other two inputs, Wg^y. = ^ 5A2^ dx - 2® • 2® = 2®. Therefore, = 24, where six 

bits are located to the right of the binary point and the two inputs from A2 are scale 
shifted six positions as depicted in figure 31. 



Figure 31.- input to A4. 

Shift pulse generator .- The shift pulse generator was developed to provide the set 
of 24 shift pulses required for operation of the accumulators. Since the accumulators 
operate at speeds of 1 /isec, the generator has to count the set of pulses and make deci- 
sions at speeds approaching the upper limit of the logic used. Carry delays associated 
with conventional counters render such a design impractical at this speed. To minimize 
delays, a feedback shift register was utilized for counting. Horton (ref. 11) has shown 
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that a 12-state counter results if a 4-stage shift register is provided with 
Feedback = AB + ABC as shown in the boxed area at the top of figure 32. This 
arrangement gives the counting sequence of table HI. 

The circuitry at the bottom of figure 32 contains the necessary controls for 
starting and stopping the generation of pulses without pulse splitting. When SW1 
is open, 12 pulses are generated as required for A1 in the input section. The 
operation of the controls is best shown in the timing chart of figure 33. The clock 
generates a 60-40 duty cycle pulse for best operation of the flip-flops. The flip- 
flops are clocked on the positive rise portion as indicated by the "action" line. 

For 12-pulse operation, the start pulse sets only M; W and the clock then set 
G at the proper portion of the clock cycle to start the generation of shift pulses. 



Figure 32.- Shift pulse generator. 
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TABLE HI.-* GENERATOR COUNTING SEQUENCE 


State 

Stage 

D 

c 

B 

A 

1 

0 

0 

o‘ 

0 

2 

1 

0 

0 

0 

3 

1 

1 

0 

0 

4 

1 

1 

1 

0 

5 

0 

1 

1 

1 

6 

0 

0 

1 

1 

7 

1 

0 

0 

1 

8 

0 

1 

0 

0 

9 

1 

0 

1 

0 

10 

0 

1 

0 

1 

11 

0 

0 

1 

0 

12 

0 

0 

0 

1 

1 

0 

0 

0 

0 


1 Me CLOCK 

ACTION 

START 

M 

W 

G 

SHIFT PULSE 
CODE 


— ^ l|isecH — 



|I« > .. L k . ^ | ^ iK« . 


( — ^ « f 



(a) 12 shift pulses. 


clock ^njinJuuLnjuuuuuuuiiTJiJiriJU^ 

START f k f 



shift pulse * — ijuiRjuTJinjiJiJi^ umiinjTJiJiJiJiJuij ( 

F = DG i \ ft 


(b) 24 shift pulses. 

Figure 33.- Shift-pulse-generator timing. 



When 12 of these pulses are generated, the shift register returns to its original pattern 
of 0000; this turns on W which immediately stops further generation. The power 
amplifier has three (positive logic) OR inputs and, hence, either W or G can inhibit 
the generation by holding the input off. 

For 24-pulse operation, SW1 is closed and, hence, the start pulse sets both D 
and M. When D is reset in the middle of the chain, it sets M which again effects 
the generation of the final 12 pulses. The gap of 2 pulses in the middle of the pulse train 
allows time for switching to take place in other portions of INTROL. The end of a 
24 -pulse operation is detected by /3p, the positive transition of the signal F = DG 
shown in figure 33(b). 

Counter 1. - The 6i \p pulses which enter A3 are integrated by Cl so that 
(Cl) = £5t// = Therefore Cl is bidirectional; it is modeled after the typical counter 
of figure 34. It is recalled that j3 = 1 and (3 = 0 signify that the angle being computed 
is angle (3 and angle ip, respectively. Therefore, the Cl operation is given by 

N[C1] = /3 © S[St//] 

where 

\ 

S[5\p] = S[A4] 



Figure 34.- Typical bidirectional counter. 
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Controls .- The controls required to effect the sequential operations of the 
ijy - computer subroutine are explained in the form of Boolean equations. Flip-flops X, 

Y, Q, and E are defined, respectively, as 

S x =|6x| R x = M X 0 F 

where M x signifies that the operation is in the X mode and F is defined in figure 33 

Sy = 1 6y | Ry = My/3p 

S Q = XY Rq = XY 

= ^X + ^Y ^E = P[^j/3pE 

The E (error), X, and Y mode signals are, respectively, defined as 

m e = e 
M X = QE 
My “ QE 

The dip increments are generated by 

M = a E + EP[A4]/3 f (22) 

The shift-pulse-generator start signal is /3j + 1 6^ | where 

% = (Y“x + Xofy^E + (X + Y)/3 e (23) 

which can be derived from the simple function I = E + XY. This arrangement takes 
care of any random occurrence of 6x and 6y increments including the simultaneous 
occurrence. The function I is also used as 

R P[A4] = % 

The mode signals are used to establish the proper operational signs and the proper 
interconnections for the four accumulators. Their logical equations are as follows: 

N[A2x] = S[6x] 

N[A2y] = Spy] 
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N[A3] = S[ 6 x]M x + S[ 6 y]M Y 


(24) 


N[A4] = S[ 6 x]M x + S[ 6 y]M Y + S[ 6 ^]M e 


(25) 


B[A3] = ((QA[A2x] k + QA[A2y] k )D[SPG2] + (qS[A2x] + QS[A2y])D[SPG2j)M E (26) 


B[A4] = ^M x A[A2y] k + M y A[A2x] k )D[SPG2] + (M x s[A2yJ 
+ M y S [A2xJ ) D[SPG2] ) 6 7 + M E AfA3] k 


(27) 


Equation (27) states that the B input of A4 is connected, during the D signal 
(i.e., during the first 12 shift pulses (fig. 33)), to the Ajj signal of A2y if the 
^-computer is operating in the X mode or to the A k signal of A2x if it is operating 
in the Y mode, and during the D signal, to the appropriate A2 sign. Before 
entering A4, the data are delayed by 6 r or 6 shift pulses. (See scale delay register 
of fig. 31.) Finally, if the i//- computer is operating in the error mode in which the 
"servo" is nulling, B[A4] is driven without delay by A k of A3. 

Binary rate multiplier .- As has been indicated, the binary rate multiplier (1) per- 
forms the Q multiplication (as required, for example, in eq. (1)) and (2) changes the 
angular unit of measure from radians to v radians (to simplify the theta decode logic). 
Since u is an irrational number and hence cannot be precisely represented in digital 
form, some precision must be established for its representation in the binary rate multi- 
plier. The resolution of \f/ was set at ~0.9° in the section pertaining to A4; this reso- 
lution represents =0.9° out of 360° or «1 part in 400. In order not to deteriorate appreci- 
ably the precision limited by this resolution, v is represented to «1 part in 1000 . 

Hence, the binary rate multiplier is designed to be a 10- stage device and follows the logi- 
cal design of figure 35 except that it is bidirectional instead of unidirectional as shown. 

It is noted in figure 35 that for I input pulses, FI pulses are emitted; 

|F| is a binary number <1 - that is, |f| = .^ 2*3 • • • f n> where fj are binary bits. 
Appropriate values of 1 * 1 , for various trajectories, are found in the binary columns of 
table IV. 
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Counters 2, 3 , and 4.- It is seen in figure 36 that |60| pulses from the binary 
rate multiplier go to the merging counter C3 if (C3) # 0; otherwise, the pulses go to 
error counter C2 which contains e = $ - 6 C . Then e is integrated in C4 which 
contains 6 = 0 C + Of. As noted in the figure any change in 6 C , whether from its own 
S/2 -1 stage or from stages ^2° of C2, is negatively reflected back to C2 to correct 
the 6 C portion of e . 

This implementation is detailed in figure 37. It is noted that C2 differs from the 
typical counter of figure 34 only in that additional inputs to each stage are provided. 

When (C2) = e < 0, the value appears in two's complement notation. Therefore, since 

6 = Se (28) 

and (C4) = 6, N[C4] = 0, always. Hence, C4 is a basic asynchronous unidirectional 
counter except for one change. The constraints associated with equation (28) require C4 
to count or add as specified in table V. Thus the values on the left of the table are such 
that e = Qc + Of, where 0 £ (6> c = 0(mod 1 octant)] < 8 and |<9f | ^0.5 octant. The right 
side of the table shows the corresponding binary number. This pattern is the counting 
sequence of a standard natural binary counter (in which /3 transitions of each stage 


MJ 



PRESET 

Figure 36.- Interconnections of C2, C3, and C4. (Counter dimensional unit is octant where 1 octant = 45°.) 
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Figure 37.- Detail of section A of figure 36. 



C2 

e 


Ck 

0 


TABLE V.- COUNTING SEQUENCE FOR C4 
p, 6 C , and Of in octants] 



trigger the next most significant stage) 
except for the 2° stage which triggers upon 
the occurrence of a transitions of the 
2"1 stage. This simple implementation 
deviation can be seen in figure 37. 

The design of C3 is very similar to 
that of C2. The addition or subtraction 
operations of the two counters are summa- 
rized in the following Boolean equations: 
For \5<p\ input, 

N[C3] = N[C2] = N[Cl] = 0 © S[A4] 

For |ae c ^| in P ut > 

N[C3] = N[C2] = S[C3] 

For the negative feedback from C4, 

N[C2] = 1 

For initialization, 

N[C3] = 1 
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Theta Resolver 


The theta resolver of figure 36 must resolve 0 C into 6x and 6y increments. 
Table VI is the theta resolver truth table. The X's in the table stand for "don't care" 
conditions. From the table it is seen that 

|6xf = ABC + ABC = BC 

therefore |6x| = B + C 

S[5X>AB + AB (bottom = X = o) 

|6y[ = ABC + ABC = BC 

therefore 1 6y | = B + C 

S[6y] = AB + AB = A 


TABLE VI. - THETA RESOLVER TRUTH TABLE 


8 C , octants 

H 

CO 

r o>’ 

A 

M 

Sjpyj 

A 

2 2 

B 

2 1 

c 

2° 

0 

0 

0 

1 

0 

0 

X 

0 

0 

1 

1 

0 

1 

0 

0 

1 

0 

0 

X 

1 

0 

0 

1 

1 

1 

i 

1 

0 

1 

0 

0 

1 

i 

0 

X 

1 

0 

1 

1 

i 

1 

1 

1 

1 

0 

0 

X 

1 

1 

1 

1 

1 

1 

0 

1 

1 


Merging Function 

An approximation to the merging function given in figure 7 is shown in figure 38. 
The abscissas are related by D = dp. The stair steps show the domains for which vari- 
ous |60 C( £| corrections are made. The graph was used to determine the values given 

in table VII from which the following Boolean equations were obtained: 


'top X = 0 \ 

^bottom X = lj 
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radians 


0.50 n radians = 2 i octants = 90 


2 octants 


50 100 


102U 

_l I 

500 1000 


D = (A5) = dp , increments 

Figure 38.- Merging function approximation. 


A " (|®«c<#.| - 2 ' 4 ) - 10 + 11 + 12 

B =(K«|,|= 2 ' 3 ) = ' 8 + 9)A 

C= (| 69 c*l = 2- 2 )-(6 + 7)AB 
D = (|O» C 0) = 2" 1 ) = « + 5) ABC 
E = (|Se c<J> j = 2°) = (2 + 3)ABCD 

F = (| 6e c<f.| - 2l ) - 155555 


+ . . . + 24 


TABLE Vn.- MERGING INCREMENTS 


Since the corrections are mutually exclusive, they 
are gated in parallel into the counters without 
carry interference. 


Most significant 

b6„ > correction 
c <p 

ON stage of A3 

Degrees Octants 

1 

90 2 1 

2, 3 

45 2° 

4, 5 

CO 

CO 

1 

H- 1 

6, 7 

lli 2 -2 

4 

8, 9 

CO 

1 

CM 

0 

1 

CO 

4 2~ 4 


Input 


The input section of S-BLOCs 2 and 3 of figure 26 contains the three natural binary 
input registers A, B, and C and the accumulator 1 which takes the various first dif- 
ferences among the three input points contained in the registers. 

Registers .- Since RA, RB, and RC are implemented in x-y pairs, six 12-stage 
standard shift registers are connected as in figure 16 so that RC shifts into RB and 
RB into RA. This "linear" shifting enables the registers to be updated by a new input 
point. The registers also operate in a "circular" mode in which the serial output feeds 
back into the input as required by Al. 

Accumulator 1 .- Accumulator 1 is represented by the typical INTROL accumulator 
of figure 28; it has 12 stages and uses the previously described shift pulse generator 
(fig. 32) which is operated with SW1 open. Since RA, RB, and RC form the various 
serial inputs to Al, their contents must be shifted to Al. To prevent destruction of the 
data during this process, they are also fed back to the registers' inputs in accordance 
with the circular mode. 

Negative quantities in the registers are represented in two's complement form 
(ref. 12) and, hence, Al operates on these quantities just as though they were positive. 
Therefore, the operation of Al for RA is always addition and for the other two reg- 
isters is always subtraction as required for the first differences of steps (9) and @) in 
figure 21. 


Data Entry 

As indicated in figure 23 two sets (12 each) of manually operated toggles are pro- 
vided for entering the x and y natural binary data. After flip-flop G has been reset 
by SW1, a main program command pulse causes the data to be transferred in parallel 
into RC. Since both set and reset inputs are connected to double-throw switches, 
resetting the register prior to entry is not necessary. In addition to the manual switch 
data entry, paper tape input is provided as an alternate mode of operation. Tape, pre- 
pared on the Flexowriter in the format of table I, is interpreted by a standard tape reader. 
Since natural binary coding is used in the input registers, the binary coded decimal (BCD) 
on the tape is converted to the natural binary code before entry into RC. Also, all nega- 
tive input points are converted to two's complement representation. 

Code detector .- The code detector is used to recognize the non- numeric or control 
characters emitted by the tape reader. In INTROL these special characters serve the 
following functions: 
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A "space" character which on tape separates the x- and y-data and controls flip- 

flop X of figure 23 which channels the x and y input data to their 
respective C registers 

(+) and (-) sign characters which control flip-flop S that "turns on" the two's com- 
plementer for negative quantities 

R carriage return which stops the tape reading and sets flip-flop X for next 

data point 

17 

Since a seven level code is utilized and only a few of the possible 2 1 combinations 
are used, many "don't care" conditions exist which permit a relatively simple synthesis 
of the detector. The particular detection codes used are specified in table I. 

BCD to binary converter .- The BIDEC method (ref. 6) is used for the BCD to binary 
conversion. As related to this particular application the BIDEC algorithm is given in the 
section "Operations of the Experimental Model." 

TABLE Vm.- TRUTH TABLE FOR "3 SUBTRACTER" 


If decade reads 

Change decade to 

Binary coded decimal 

Equivalent 

decimal 

number 

Equivalent 

decimal 

number 

Binary coded decimal 

D 

2 3 

C 

22 

B 

2 1 

A 

2° 

D 

2 3 

c 

2 2 

B 

2 1 

A 

2° 





0 to 7 

(No change) 





1 

0 

0 

0 

8 

5 

0 

1 

0 

1 

1 

0 

0 

1 

9 

6 

0 

1 

1 

0 

1 

0 

1 

0 

10 

7 

0 

1 

1 

1 

1 

0 

1 

1 

11 

8 

1 

0 

0 

0 

1 

1 

0 

0 

12 

9 

1 

0 

0 

1 

1 

1 

0 

1 

13 

10 

1 

0 

1 

0 

1 

1 

1 

0 

14 

11 

X 

X 

X 

X 

1 

1 

1 

1 

15 

12 

X 

X 

X 

X 


The specialized "3 subtracter” is synthesized by finding the logic which implements 
truth table VIII. Since counts 14 and 15 never occur, they are considered "don't care" 
conditions. From the reduction maps of figure 39, the following implementation for the 
four flip-flops of each decade is obtained: 


65 










s D = 0 

R D = AC + BC = c(a + b) 

Sq = DC (a + b) = DRj) = /3 q 
Rq = CD 


S B = adb 

R b = ADB 


Hence, 

T B = AD 

s A = ad 

R A = AD 

Hence, 

t a = d 

It is noted that the characters are emitted from the tape reader with the most sig- 
nificant character leading and that BIDEC requires the characters to be shifted by RT 
of figure 23 with the least significant character leading. Since the shift register used is 
unidirectional, the required rearranging is accomplished as follows: 

(1) Connect (with fictitious switch L) the serial output of the shift register to its 
serial input. 

(2) Shift register 12 times. 

(3) Read the BCD character, 4 bits in parallel, into the units position as shown in 
figure 23. 

(4) Repeat steps (2) and (3) for a total of four characters. 

Since RT is 16 bits long, its contents after reading the four characters are 
arranged as shown in figure 23. This arrangement is now ready for the BIDEC operation. 

Two's complementer .- If the absolute value of a negative number is subtracted from 
zero in the typical INTROL accumulator of figure 28, the two's complement results. 

Thus, since the two's complementer may be considered an accumulator in which A k = 0 
always, equations (17), (20), and (21) reduce, respectively, to 

A k+ 1 = B k P k + B k P k = B k © P k (29) 
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Sp = NBk 


(30) 


Rp = NB k = 0 (31) 

For the complementation of each number, P is initially reset. Equations (30) 
and (31) state that for complementing (i.e., N = 1) P is set the first time that B k is 
true in the serial operation and thereafter P remains set. Upon this hypothesis, equa- 
tion (29) states that for complementing, A k+ j = B k © 0 = B k until the first B k = 1, 
after which A k+1 = B k © 1 = B k . For positive numbers (i.e., N = 0), P never gets set 
and, hence, A k+ j = B k . Therefore, the B k number does not get changed. 

EXPERIMENTAL RESULTS 

In general, the experimental results are in good agreement with the theory and 
confirm that the mathematical technique and implementation thereof are feasible. The 
curves presented in this section are reproductions of those actually generated by the 
model and were selected because they emphasize both merits and weaknesses of the 
present model. 

Evaluation of Generated Curves 

The coordinates of the input data points of figure 1 were determined and read into 
the experimental model. The generated curve of figure 40 shows close agreement with 
that of figure 1. 


Y 



0 100 200 300 400 


Figure 40.- Points of figure 1 interpolated by experimental model. 
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Figure 41.- The figure eight. 

A considerably more difficult task is exhibited in figure 41. The figure- eight curve 
demonstrates the ability to 

(1) draw smooth curves through widely spaced points 

(2) operate on nonequal argument data 

(3) produce multivalued functions 

(4) process data in all four quadrants 

(5) handle infinite tangents 

(6) produce both positive and negative rotations 

At first glance, the curve between points 8 and 9 appears to be in error. However, when 
the point sequence indicated in the figure is carefully considered, the distorted shape is 
found to be a reasonable interpolation. 

The. curve between points 8 and 9, however, does reveal a weakness of the experi- 
mental model. The abruptness in the turning is a result of the inability of its merging 
hardware to handle effectively any points which are spaced >0.7 inch. This range, how- 
ever, is easily extended by adding similar stages of logic. 


I 
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Figure 42.- Merging examples. 


Y 



Figure 43.- Circle through two points and one slope. 


Figure 42 shows some results of merging 
tests designed to insure satisfactory merging in 
all directions. The point separation of ~0.7 inch 
was chosen since, as previously mentioned, this 
distance represents the maximum effective merging 
range for the model. 

As predicted in the section "Derivation of 
Output Function," the circle of figure 43 does have 
four sides slightly flattened. It is recalled that 
this tendency to flatten results from a mathematical 
approximation and that for normal curve interpola- 
tion this flatness is insignificant. This particular 
circle was drawn clockwise, first from A to B 
and then from B to A. In addition to the 
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coordinates of A and B, the input data included the slope at B for the A to B 
portion and the slope at A for the B to A portion. 

The curves in figures 40 to 43 were generated in the circular mode. The cardioid 
mode is demonstrated in figure 44. The curve was generated first from A to B over 
the lower path and then from A to B over the upper path. The input data consisted 
of the coordinates of A with respect to a coordinate system centered at B and alined 
with the slope (zero) of B. The upper or lower path is selected by integrating xfy over 
a path in the upper or lower half-plane, respectively. An explanation of this selection 
can be seen in basic equation (1), 

0’ = Q^' 


by observing the initial <p' for both ip = +ir and \p = -it. Furthermore, since 
d <p' = Q dt//' and Q > 0, the sense of rotation of is the same as that of t//' . In 
general, then, as a point moves along the interpolated path, |t//| and |<£'| = Q|t//'| 
decrease monotonically, both ending with zero slope. 

Careful inspection of this cardioid curve and the circular curve reveals a slight 
discrepancy in symmetry. This distortion results primarily from the lack of smoothing 
(ref. 6) in the binary rate multiplier. For applications in which it is important to correct 
this symmetry, as perhaps in the control of a machine tool, two or three stages of 
smoothing are recommended. A secondary contributor to the distortion may be the 
coarse resolution in several components of the system where, in the interest of simple 



Figure 44.- Cardioid through two points and one slope. 
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and minimal hardware, the model was designed with fewest stages. For example, some 
of the computations are resolved to only ~3°. In some applications, however, this resolu- 
tion should be extended by a factor §4. 

Finally, the lemniscate mode is demonstrated by the curve of figure 45. Again two 
points and one slope dictated the path which was drawn. It is noted that the entire path 
lies in the same quadrant with the starting point A. The Y-axis then represents an 
unstable line on which the lemniscate diverges. Thus |^'| must be limited to ir/2. 


Y 



Figure 45.- Lemniscate through two points and one slope. 


i//' Boundaries 

In general, the ijs' limit occurs when reaches a value such that 

| <j>' - xp'\ = 7 r. This limit may be stated as a function of Q as follows: 

it> \4>' - *f/'\ = |Q<K - ^'| = |^'(Q - 1)| 

For Q > 1, 

|*'(Q- l)| = M(Q- 1) 

and, therefore, 

M'fQTi) < 32 > 

To cover all four quadrants, xj/' must range over 27r. This range can be effected 
with a maximum value of |i//' | of ir. Substituting this maximum value in equation (32) 
gives 
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Hence, 

(Q - 1) < 1 

and, therefore, 

Q < 2 

Thus, any two points in the finite coordinate system can be connected and termi- 
nated with any prescribed slope if 1 < Q < 2. Note that the interval is open. Therefore, 
for the important circular mode in which Q = 2, an unstable though somewhat trivial 
domain exists: the line for which \p' = ±ir. 

Computer Error 

It is recalled that the design of the ^/-computer is based on the approximation 
6\f/ = tan 6i p. Near the origin, however, 6\p can be as great as 45° per increment; this 
45° angle would be computed as *57°. This discrepancy of 12° is prohibitive. In the 
model, it was reduced to an error of <1° simply by presetting accumulator A3 to O.OI 2 
before the computation of each angle. 

Initialization of (C3) 

In step (D of the operations algorithm (p. 32), each new segment l of interpola- 
tion is initialized to -0 C , the angle at which the last step entered a data point. In general, 
this angle is not _ j, the "average" slope of the (l - l)th segment into the point. Since 
the initialization should be and not 0 C , prohibitive angular discontinuities result. 

Therefore, for the curves in figures 40 and 41 each initialization was manually adjusted 
to when |^c“^-l| was large* 

The correct initialization may be implemented by using Cl which contains, at the 
end of a segment, a measure of the desired angle (3^ _ j. Unfortunately the unit of (Cl) 
is radians and, hence, the quantity must be counted into the binary rate multiplier for con- 
version to ts radians and then recounted into C3. Step (2) of the operations algorithm is 
corrected by the following procedure: 

(I) Reset C2 and C3. 

(D Set (C2) = -(C4 c ). 

(3) Count-down (Cl) through the binary rate multiplier into both C2 and C3 
and apply appropriate multiplying factor and operational signs to leave (C2) = j - 0c 
and (C3) = -/3j _j. 

(?) Integrate (p and simultaneously count the angle into C3 to leave 
(C3) = <p - j3 z _ r 

(5) Begin the interpolated segment. 



RESEARCH BYPRODUCTS 


During the development of the subject controller, several innovations or research 
byproducts ensued. These byproducts are essentially separate entities which should find 
application in other areas, devices, or systems. The most important byproducts are as 
follows: 

^-computer: A special purpose incremental digital computer for integrating the arc 
tangent angle as a function of various finite 6x and 6y input increments 
The described design can be easily modified to obtain virtually any accuracy and capacity. 
If the ^-computer is utilized as a finite difference device, its speed is great. 

Bidirectional binary rate multiplier: A device which eliminates the "hysteresis" inher- 
ent in conventional binary rate multipliers when utilized in bidirectional applications 

Adder with round off: A digital accumulator which outputs rounded-off information to one 
precision but retains a finer precision for subsequent calculations 

Although developed in a natural binary system, the scheme is applicable to other weighted 
codes provided only that the weight of the least significant bit of the rounded portion is 
twice that of the adjacent lower order bit which, on turning on, carries into it. For 
example, in a BCD code of weight 1-2-4- 5, the least significant bit (bit 1) of one decade 
has twice the weight as that of the most significant bit (bit 5) of the adjacent lower order 
decade; hence such a code is suitable. 

High speed shift pulse generator: A device which, on command, generates a fixed number 
of fully formed pulses at a rate equal to the nominal speed of the logic with which it is 
implemented 

Although developed with the fixed number as 24, the scheme is adaptable to virtually any 
quantity. 

Serial two's complementer: A simple digital device for serially converting a sign- 
magnitude negative number to its two’s complement representation (and vice versa) 

NEED FOR ADDITIONAL RESEARCH 

The experimental results indicate that the merging- function — constant-Q concept 
developed in this research works well. Nevertheless, the concept was an expediency 
which does not fit harmoniously into the system. Therefore, the scheme which follows 
seeks to satisfy this need. 

From the fundamental relation (p ' = Qi p 1 (eq. (1)), 


<t> - 0 = Q(<// - 0) 
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At the beginning of a segment, 


4> = a 


and 


* = M A,B) 


Therefore, 

Q = Q(A,B ,a,$) 

that is, Q may be considered a function of the four independent variables A, B, a, 
and j3. This characteristic is demonstrated in figure 46 where, for the specific points 
A and B, a Q-family of curves is drawn for various values of ol. If the coordinate 
system of this figure is considered to be the transformed system, then the /3 = 0 slope 
is generalized to any slope by the transformation itself. Thus, for any two points A 
and B with slopes a and j8, respectively, a value of Q is found for use in basic 
equation (1) to enable generation of a smooth segment. The catch is that the smooth 
curve is not always the desired one. In general, the desired curve results if the direc- 
tion of a is such that ji//’ | decreases. 

This condition is satisfied in substantially 
all segments usually encountered. For an 
S- shape segment, however, the condition 
fails and the generated segment is not the 
desired connection. 

One approach is to ignore the failure 
condition since it seldom occurs. Another 
approach is to restrict the input data spacing 
to prevent the failure condition. A third and 
best approach is to find a simply implement- 
able solution which copes effectively with the 

failure condition and gives the desired path. Figure 46 ._ Q . fami | y connections from a to b. 



CONCLUDING REMARKS 

The prime objective of this research was to evolve a simple controller for auto- 
matically generating a smooth curve through an ordered set of data points even for 
unequal arguments, closed contours, and raw data. The "figure of merit" for the gen- 
erated curve was to be based on how well it resembled one properly drawn with french 
curves. Tacitly assumed was that sufficient input data were given to represent faithfully 
the curve without resorting to contextual information, such as characteristics of the data 
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source (e.g., frequency content, bandwidths, or intentional discontinuities), or recognizing 
that the input data fall into some analytic function or geometric shape. The effect of 
context was forcibly demonstrated when the figure eight was redrawn in the context of two 
circles; the resulting curve was significantly less distorted. Therefore, for the figure 
eight the sufficient input data assumption was violated. The foregoing and other experi- 
mental curves do indicate that given sufficient input data, the controller generates curves 
of equal fidelity to the skilled draftsman using french curves. Another way of stating 
this fact is: if enough input data are given so that the same curve results when drawn 
(out of context) by several skilled draftsmen, then the controller will also produce the 
same curve. The word "same" does not require a quantitative definition because even a 
highly subjective interpretation is perfectly adequate if consistently applied to both 
draftsmen and controller. 

Considering the magnitude of the task which it performs, the controller satisfies 
also the secondary objectives of simplicity and speed. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 6, 1966. 
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APPENDIX 


DERIVATION OF STEERING ANGLE 4> 

In figure 1 is a set of input data points a, b, c, d, e through which a smooth 
curve is drawn. Assume that the first portion of the curve is transformed to the new 
(primed) coordinate system of figure 2 such that the new origin is centered at point b 
and that the positive X f -axis is alined with /3, the slope of the curve at b. This special 
case is analyzed in the more convenient transformed system and then generalized by 
means of the transformation itself. This procedure is particularly attractive since most 
of the parameters of the problem are angles which are invariant under translations and 
are simply rotated under coordinate rotation. These properties can be seen in figure 47 
where under translation (from system 1 to 
system 2) the line pb does not rotate and the 
X-axis reference direction does not change; 
hence is invariant. However, under 
coordinate rotation (from system 2 to sys- 
tem 3), the reference direction changes by 
the amount of the rotation; hence i// changes 
by an equal amount and, thus 

= * - 0 (33) 

Another useful property is the invari- 
ance under rotation (and of course under 
translation) of a change in t//. Thus, since 
/3 is constant for x a ^ Xp < x^, 

di// 1 = dt// - d/3 = d^ (34) 

Consider the present problem to be that of determining the angle associated 
with the tangent of figure 2. The angle can then be transformed to <fi which can 
be used to guide point p from a to b. 

Since b f is a destination of p T , the direction of b T with respect to p r is 
apparently an important parameter. The angle i^ T is therefore defined, as indicated in 
figure 2, as the instantaneous angle of the line connecting b T and p T as p’ moves 
along the curve — that is, 

- * mr %) 

where x T and y T are the coordinates of p r . 



Figure 47.- Effect of transformation on ift. 
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Let 0' = 0'(0') 
(eq. (D) 


and, in the interest of obtaining simple relationships, assume that 
0 ’ = 00 ' 


where Q is a constant. Several significant cases then result. 

Case I — Q = 1: Obviously case I is the linear case in which 

0 ' = 0 ' = tan -1 ^!^ = c 

where c is a constant for x^ s X ' < x^. Hence 0' would guide p' over a constant 
slope or linear path from a' to b*. 

Case II - Q = 2: From equation (1), 


0 ’ = 20 ' 

In figure 2 it is seen that ^ = tan 0'. Therefore, 

= tan 20' = -3- tan = Kx) 
dx 1 - tan 2 0' x _ m 2 



(35) 


a first order, homogeneous differential equation. (Prime with x and y is implicit.) 
Let 


then y = vx and 


Separating the variables gives 


v 


y 


x 


dy dv 

dx dx 


2v 

1 - v 2 


x 

dx 


2v 


1 _ v 2 


- v 


dv 


2v - v + v 3 
dv(l - v) 2 


Therefore, 

dx _ (l - v 2 )dv 
x = v(v 2 +l) 

By use of partial fractions, 

1 - v 2 _ A B 
v(v 2 + l) v V 2 + 1 


(36) 


(37) 


(38) 
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and therefore, 


A = 


+ 1 


v=0 




Substituting A = 1 


Hence, 


in equation (38) gives 


1 - _ 1^ B _ v 2 + 1 + Bv 

v(v 2 + l) v v 2 + 1 v(v2 + l) 


B = -2v 


Substituting the values of A 
tion (37) results in 


Integrating gives 


and B into equation (38) and equation (38) into equa- 



ln x = In v - In v 2 + 1 + In c 


and, hence, 



i) 


= In c 


Taking the antilogarithm gives 

±c «x(± + v ) 

Therefore, from equation (36), 



By substituting 2r for c, 

x 2 + y2 ± 2yr = 0 


and thus 

x 2 + (y ± r) 2 = r 2 


Therefore, the resulting path for <p ' = 2i[/' is one of 
a family of circles or radii r, centered on the Y-axis and 
displaced from the origin in both positive and negative direc- 
tions by a distance equal to their radii. (See fig. 48.) 


Y' 


X' 



Figure 48.- Circular paths for 
4 )’ = 2 i/i' . 
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Case HI — Q = 3/2: From equation (1), 

< p ' =!^ t 

Polar coordinates will be used to show 
that case in also yields a familiar geometric 
curve. First, however, it is convenient to 
introduce the lemma 

tan 6 = r ~ 
dr 

where r and 9 are the standard polar 
variables and 6 is the angle between the 
radius vector and the tangent line (fig. 49). 


Y* 



Figure 49.- Relation of 6 to O' , ip' , and 8. 


For the proof of the lemma, assume that the curve is given in the form r = r (9), 
where r(0) is differentiable. From figure 49, it is seen that 

tan 6 = tan (<£' - x[/') 

= tan(0' - 9 + 7r) 


Therefore, 


tan 6 = ta n(4>' - 9) 

_ tan <ft* - tan 9 
1 + tan 4 > ' tan 9 


dy _ y 

_ dx x 

_ x dy - y dx 
x dx + y dy 


Recall from the polar to Cartesian transformation that 


and, therefore, 


x = r cos 9 
y = r sin 9 


dx = -r sin 9 d0 + cos 9 dr 


(40) 


(41a) 

(41b) 


(41c) 


dy = r cos 9 d9 + sin 9 dr 


(4 Id) 
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Substituting equations ( 41 ) into equation ( 40 ) gives 

q _ r cos 0(r cos 0 d0 + sin 0 dr) - r sin 0(-r sin 0 d0 + cos 0 dr) 
r cos 0(7 r sin 0 d0 + cos 0 dr) + r sin 0(r cos 0 d0 + sin 0 dr) 

r^cos^0 d0 + r^sin^0 d0 r^d0 d0 


which proves the lemma. From this lemma and figure 49 , 


r ™ = tan 6 = tan(i^'/2) = tan ^ = -cot(0/2) 
Separating the variables results in 


dr _ d0 _ _ sin(0/2)d0 
r cot(0/2) cos(0/2) 


Since the numerator is the differential of the 
denominator (except for a factor of 1/2 on 
the right-hand side) equation ( 42 ) is easily 
integrated to give 

In r = 2 In cos(0/2) + In 2a 

where r > 0, cos(0/2) >0, a > 0, and 

In 2 a is the integration constant. Taking 
the antilogarithm gives 

r = 2a cos^(0/2) (- w < 6 < ir) 

and, therefore, 


r = a(l + cos 0) ( 43 ) 

the classical form of the cardioid of 
figure 50 . 

Case IV — Q = 3 : From equation ( 1 ), 


Y' 



Figure 50.- 


Cardioid for O' = 


= 3i y 

Again, from the lemma and figure 49 , 


r ^ = tan 6 = tan(2(0 - it )) = tan 20 



APPENDIX 


Y' 


Separating the variables results in 

dr _ dfl _ (cos 26>)d0 
r tan 20 sin 26 V ’ 

which is easily integrated to give 

In r = i ln(sin 26) + In a 

where r > 0, sin 26 > 0, a > 0 (fig. 51), 
and In a is the integration constant. Taking 
the antilogarithm yields 

r^ = a^sin 26 (45) 



the classical form of the lemniscate (ref. 13) Figure 51.- Lemniscate for <t>' = 3 ip'. 

of figure 51. 


Since most of the equations of the foregoing figures are recognizable more easily „ 
in polar coordinates, a polar generalization, which includes cases II to IV (case I is con- 
sidered to be trivial), is now derived. The relations 


and 

from figure 49 with equation (1) 


6ef-f 
l//' = 0 - IT 
= Q^’ 


give 


6 = (Q - 1)(0 - it) 

Therefore the lemma can be expressed as 


r || = tan[(Q - 1)(0 - ir) + irk] 


where k is any integer. Separating the variables results in 

d r _ cos[(Q - 1)(6> - it) + irk]d0 
r sin[(Q - 1)(0 - ir) + 7rk] 


which is easily integrated to give 

In r = 1 In sin[(Q - 1 )(6 - ir) + 7rk] + — In c 

Q " 1 Q - 1 

where r > 0, sin[(Q - 1)(0 - ir) + irkj > 0, c > 0, and — - In c is the integration 

Q - 1 

constant. Taking the antilogarithm yields the polar generalization 
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rQ _1 = c sin[(Q - 1)(0 - v) + 7rkj ^1 - < © < it ^1 - (46) 

where the domains of 0 result from the conditions that 

^sin[(Q - 1)(0 - 7r) + Trk] > o) =} ^0 + 2irm < [(Q - 1)(0 - tt) + irkj < it + 2irmj 

for any integer m, and 27 rm c kw. 

For case II, in which Q = 2, equation (46) becomes for k = 1 

r = c sin 0 (0 < 0 < ir) 

and for k = 0 

r = -c sin 0 (ir < 0 < 2ir) 

which agrees with equation (39), as can be seen from the following transformation to 
Cartesian coordinates: 

_ r _ r 2 x 2 + 

~ sin 0 ” r sin 0 y 

For case HI, in which Q = 3/2, equation (46) becomes for k = 1 

r 1 /2 = csin g + |j = C cos(|j 

Squaring both sides gives 

r = c 2 cos 2 ^ = c 2 ^(1 + cos 0) 

2 

which agrees with equation (43), where ^- = a. 

£k 

For case IV, in which Q = 3, equation (46) becomes for k = 2 and k = 0, 
respectively, 


(-ir < 9 < ir) 


r 2 = c sin 20 


o<0<!> 

‘* <9< f 


which agrees with equation (45), where c = a 2 . Also for k = 1 and k = -1, 
respectively, 

/ | < 0 < 7T 

r 2 = -c sin 20 « 

\¥< e<2 ' 
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Thus, equation (46) becomes 

= c sin [(Q - 1)0] (47) 

for integer values of Q > 0 and for k = Q - 1, and 

r 2(Q-l) =c 2 sin 2[ (Q . 1)e + |j 

= C 2 COS 2 [(Q - 1)0] 

= c’ |l + cos[2(Q - 1)0]| (48) 

c 2 

for integer values of (Q + 0.5) > 1 and for k = Q - 0.5 and c' = 

Other popular geometric shapes can probably be found by solving equation (47) for 
additional integer values of Q or solving equation (48) for additional "half-integer" 
values. Also, Q may be quartered or further subdivided to obtain other shapes. With- 
out attempting to exhaust all of the important shapes, however, a significant conclusion 
is stated in the following theorem which is based on figure 49: 

Theorem: If 0 is the polar angle of a point on a curve and cp is the tangent 
angle of the curve at the point, then for -ir = (\p = 6 - it) < ir simple linear relationships 
exist between cf> and rp, such that p generates smooth well-known geometric shapes 
which go through the coordinate origin at zero slope. 

All that remains to be proved is the origin zero- slope portion. That the origin is 
a solution of equation (46) is easily seen by substitution; its slope may be determined as 
follows: 

lim rQ"l = 0 = lim c sin[(Q - 1)(0 - it) + irk] 
p— 0 p—0 

where p — 0 means "as point p approaches the origin." Therefore, 

lim Kq - 1)(0 - it) + irkl = irm 
P—0 L 

where m is any integer. Letting m = k gives 

lim j(Q - 1)(0 - tt)1 = 0 
P—0 

and since (Q - 1) * 0 

lim 8 = it 
P—0 

Therefore, the curve goes through the origin at zero slope. 
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